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Abstract. It has been shown that a discontinuity in the derivatives of the sound speed at the edge of the convective 
regions inside a star gives rise to a characteristic oscillatory signal in the frequencies of stellar oscillations. This 
oscillatory signal has been suggested as a means to study the base of the outer convection zone in low mass stars 
and possibly the outer edge of the convective core in high mass stars. Using stellar models we show that because 
of a phenomenon similar to aliasing in Fourier transform, it may not be possible to use this signal to detect the 
convective core. Nevertheless, it may be possible to determine the size of convective cores using the frequency 
separation fn+i^e — fn/. 
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It has been shown that a discontinuity in the deriva- 
tives of the sound speed in the stellar interior, like those 
arising at the edge of a convective region, introduces a 
characteristic oscillatory signature in the frequencies of 
low degree modes (Sough 199C|; Monteiro, Christensen- 



Dalsgaard & Thompson 1994; Roxburgh & Vorontsov 



% 



1994; Basu, Antia & Narasimha 1994). From this os 



dilatory signature it has been possible to put limits 
on the extent of overshoot below the solar convection 



zone (Basu 1997 ). Monteiro, Christensen-Dalsgaard fc 



Thompson (200C ) have pointed out that this oscillatory 
signature can be used to study the location of the base 
of the convection zone as well as the extent of overshoot 
below this base in other stars using asteroseismic data for 



only low degree modes. Monteiro, Christensen-Dalsgaard 



fc Thompson (1998 ) have also suggested that the same os- 
cillatory signal can be used to study the size of convective 
cores in massive stars. However, it is not clear if such a 
signal has been demonstrated using actual stellar models 
with convective cores. 

Similarly, helioseismic inversions for the rotation rate 



in the solar interior (Thompson et al. 1996; Schou et al 



199S) have shown the presence of a tachocline (3piegel & 
Zalm 199^) around the base of the convection zone, where 



the rotation rate varies from differential rotation inside 
the convection zone to almost solid body like rotation in 
the radiative interior. This sharp transition in the rotation 
rate introduces an o scillatory signature in th e frequency 
splitting coefficients ( Mazumdar fc Antia 2001 ), which can 



be used to study the characteristics of tachoclines in other 
stars. 

In this work we attempt to study the possible oscilla- 
tory signal arising from convective cores in stars more mas- 
sive than the Sun. However, we do not find any significant 
oscillatory signature of the core and try to identify the rea- 
son for the apparent absence of this signal. We also study 
how the amplitude of oscillatory signal from the outer con- 
vection zone varies with stellar mass and age. Since even 
the signal from the outer convection zone is rather weak, 
such studies may help us in identifying promising stars 
where attempts can be made to detect such an oscillatory 
signal. In the absence of a substantial oscillatory signal 
from the convective core of stars, we try to identify alter- 
nate means to detect the presence of a convective core and 
to measure its size. 

2. The technique 



In the asymptotic limit (Christensen-Dalsgaard fc 



Berthomieu 1991) the eigenf unctions of stellar oscillation 
modes have the general form cos{lut + 0), where r is the 
acoustic depth given by 

T — I -——- dr 



1 

c{r) 



(1) 



Here R is the stellar radiu s. Using this form for the eigen- 
functions it can be shown ( Monteiro et al. 1994 ) that if we 
have a discontinuity in the derivatives of the sound speed 
at a point corresponding to an acoustic depth Tc then the 
frequencies of p-modes can be written as 
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n,l = ^n.e + M^nj) sin(2w„_fTc -f ( 



(2) 
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where n is the radial order and i the degree of the mode 

(s) 

and a;„ ^ is the smooth part of the frequency arising from 
the smooth variation in the sound speed with depth, while 
the second term is the oscillatory component. The ampli- 
tude A(uj) is a smooth function of w. Here we have ne- 
glected terms involving the degree £ since their effects are 
small for low degree (£ < 3) modes. In this work we are 
primarily interested in stars other than the Sun, where it 
is not possible to detect oscillations with higher degree. 
Similar oscillatory terms can arise in the splitting coeffi- 
cients if there is any steep variation in the rotation rate 



such as what happens in the tachocHne region ( Mazumdar 
& Antija 200 1|). 



Following Basu et al. (1994 ) , we take the fourth differ- 
ence of the frequencies with respect to n to enhance the 
oscillatory signal. Another advantage of taking the fourth 
difference is that the smooth part of the frequencies be- 
comes negligible and we do not need to include it in our 
analysis. This will of course, depend on the smooth com- 
ponent of variation of the sound speed, but at least for 
the Sun this component is found to be negligible. The 
sharp variation in sound speed inside the second heUum 
ionisation zone also gives rise to an oscillatory signal, but 
in this case the region with reduced adiabatic index, Fi, 
has a small but finite width and in some sense it can be 
considered as two steps of opposite sign located close to 



each other. Monteiro & Thompson (1998) have shown that 
the amplitude of such a signal would be modulated with 
an oscillatory factor of the form sin^ {^P) , where /3 is the 
acoustic half- width of the Hen ionisation zone. We have 
not included this factor in our adopted form for the oscil- 
latory signal. Thus for the Sun there are two oscillatory 
components - one arising from the base of the convec- 
tion zone and the other from the second helium ionisation 
zone. A similar form may be expected for other stars with 
M < 1.2 A/©, where there is no convection in the core. 
For more massive stars with convective cores we may ex- 
pect a third component arising from the boundary of the 
convective core. 

Since there are no observations yet of stellar oscilla- 
tions with sufficient accuracy to detect this oscillatory sig- 
nal, in this work we construct stellar models with different 
mass and age and calculate their oscillation frequencies. 
Using the computed frequencies we attempt to check if 
the location of these regions of rapid variation in sound 
speed can be determined using the frequencies of only low 
degree modes, which are likely to be observed by forth- 
coming asteroseismic missions. 

Apart from frequencies, we also study the oscillatory 
signal in the splitting coefficients arising from possible 
tachoclines in stars. There is no theory to determine the 
position of tachoclines in stars, but for simplicity we gen- 
erally assume that the tachocline is located at the base 
of the out er convection zone, w hich is in fact the case for 
the Sun ( Spiegel fc Zahn 1992| ). However, in some cases 
we have studied the effect of varying the position of the 
tachocline to see how the oscillatory signal varies with its 



depth. The main advantage of using the splitting coeffi- 
cients arising from an assumed tachocline is that, in this 
case there is only one layer where the rotation rate changes 
rapidly and its location can be altered to study how the 
signal varies. This exercise is not possible with frequencies 
where there are multiple layers with rapid variation in the 
sound speed and their locations cannot be easily shifted 
over the entire stellar radius. 

In order to study the oscillatory signature in splitting 
coefficients we assume a model tachocline rotation profile 



of the form (Mazumdar & Antia 2001) 



n{r,e) 



(517(5 cos^6>-l) 
1 -I- exp[(rd - r)/w] ' 



(3) 



where SQ is the extent of variation in the rotation rate 
across the tachocline, lying at a radial position of rj, and 
having a half-width of w. Using this model rotation profile 
we can calculate the corresponding splitting coefficients 
for a star. Because of the choice of latitudinal dependence 
only the splitting coefficient C3 (n, i) is found to be non- 
zero. Neglecting the smooth part, the fourth difference can 
be written as 



6^C3inJ)^ Uo + ^ + 






OJ, 



0.2 
2 



smi2uJn^iT - 



(4) 



The parameters ao,ai,a2,T and (j) can be determined by 
a nonlinear least squares fit. For the frequencies also we 
use the same form for fitting except that there are mul- 
tiple oscillatory terms. Each oscillatory part is defined by 
a similar set of parameters, but the number of parame- 
ters to be fitted is larger and hence it is more difficult to 
obtain a proper fit. For distant stars it will be possible 
to detect modes with £ = 0, 1, 2, 3 only and hence we use 
only these modes in our study. As mentioned earlier, the 
amplitude of the signal arising from the He 11 ionisation 
zone has an oscillatory factor and may not be properly 
represented by the form we have adopted. However, the 
half-width /3 being quite small, the 'wavelength' of oscilla- 
tion in amplitude would be quite large and in the limited 
frequency range that is likely to be observed we expect 
less than half a wavelength to be accommodated. Such 
a variation can be approximated reasonably well by the 
three terms that we have incorporated. Hence we have lit- 
tle difficulty in fitting the oscillatory signal from the He 11 
ionisation zone. 

To illustrate the oscillatory signal. Fig. n] shows a fit 
to the fourth difference of frequencies for a star with 
M = 1 Mq and age 4.55 x 10^ years. This model is slightly 
different from the current solar models, since diffusion of 
helium and heavy elements in the radiative interior is ne- 
glected. One can see the two oscillatory components in 
this figure. The slowly varying component with t w 676 s 
arises from the second helium ionisation zone and the 
value of r is approximately equal to the acoustic depth 
of this ionisation layer. The more rapidly oscillating com- 
ponent (r w 2176 s) arises from the base of the convec- 
tion zone. The amplitudes of all oscillatory components 
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(Iglesias & Rogers 1996) opacity tables; nuclear reaction 
rates from the NACRE compilation (Angulo, Arnould & 



Rayet 1999| ; [Morel 1999D ; mixing length of l.SiJp, where 
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!/ (mHz) 



4.0 



Fig. 1. Fourth difference of the frequency i/, for a solar 
mass star at an age of 4.55 Gyr. A fit to the data using 
a function having two oscillatory components of the form 
described in Eq. (Q) is shown. Each oscillatory component 
is also shown separately in the hgure. 



Hp is the pressure scale height. 

We have considered stellar models with M /Mq — 
0.5, 0.75, 1, 1.25, 1.5, 1.75, 2, 2.5, 3, 4 and 5 at different ages 
spanning the main sequence phase of evolution. The 
ZAMS chemical composition of all the models were taken 
to be solar-like (X = 0.708, Z = 0.019). The characteris- 
tics of some of the representative models are summarised 
in Table ll|. Here, X^ is the hydrogen abundance at the cen- 
tre. Vac — c/ {At: Hp) is the acoustic cutoff frequency at the 
surface, tq is the acoustic radius of the star or the sound 
travel time from the centre to the surface. The radial dis- 
tance as well as the acoustic depth of the outer edge of 
the convective core, base of the convective envelope and 
the Hen ionisation zone are also given in the table. 

For each of these models we calculate the frequencies 
of p-modes for (, = 0,1,2,3 and use them to study the 
oscillatory signal using the second or fourth difference. 
The oscillatory signal also depends on the extent of over- 
shoot below the convection zone (Montciro et al. 1994) as 
well as the treatment for calculating the diffusion of he- 
lium and heavy elements in the radiative regions (BasL 



& Antia 1994). From the solar studies it has been known 



that the overshoot below the base of the convection zone is 



decrease with frequency. At still higher frequencies where 
the modes are not trapped below the solar surface the ef- 
fective value of T for the fits may also change. We do not 
include such modes as they are not likely to be observed 
and in any case their interpretation is difficult. Thus for all 



quite small (Basu 1997) and further the amplitude of the 
observed signal is consistent with that in a solar model 
without overshoot and with very little composition gra- 
dient at the base of the convection zone. The absence of 
composition gradient is probably due to some mixing pro- 



stellar models we generally restrict ourselves to frequencies 
less than the acoustic cutoff frequency tOac ~ 2B^ wherFT" 
is the sound speed and Hp is the pressure scale height near 
the stellar surface. At low frequencies the fourth difference 
generally shows steep variations, probably because the low 
radial order modes do not conform to the asymptotic rela- 
tions used in obtaining the oscillatory form. Thus we ne- 
glect these modes too while studying the oscillatory signal. 
The frequency range studied for different stars varies with 
mass and age and in order to compare the respective am- 
plitudes of the oscillatory signal we scale the frequencies 
by y^A'I/R-^, where M is the stellar mass and R the radius 
in terms of solar values. Thus we compare the amplitudes 
of oscillatory signal at frequency vq \/p/pq, where vq is a 
fixed value chosen for the Sun and p is the mean density. 



3. The stellar models 

In order to study the variation in the oscillatory signal 
for different stars we construct models for stars with dif- 
ferent masses a t various ag es using the stellar evolution 
code CESAM (JMorel 1997|) . For simplicity we have ne- 
glected diffusion of helium and heavy elements in stel 
lar interior. We use the EFF ( Eggleton, Faulkner fc 



cess operating in that region (Richard et al. 1996; Brun 



Turck-Chieze, fc Zahn 1999 ). Thus in this work we ne 



Flannery 1973) equation of state with Coul omb correc- 
tion ( Christensen-Dalsgaard fc Dappen 1992|) ; the OPAL 



gleet both diffusion and overshoot. Of course, it is possible 
that in some stars the overshoot may be larger, but that 
will only increase the amplitude of the oscillatory signal. 
Similarly, a gradient in helium abundance below the outer 
convection zone also tends to enhance the oscillatory sig- 
nal. Thus by neglecting overshoot and diffusion we will 
get a lower limit on the amplitude of the oscillatory sig- 
nal. For massive stars where the outer convection zone 
becomes shallow, it is quite likely that extent of overshoot 
as well as composition gradient may be larger than that 
in the Sun, because the efhciency of convection would be 
lower ( ^aikiaet al. 2000| ). 

Apart from this, for most of these models we also calcu- 
late the splitting coefficients assuming a model tachocline 
profile given by Eq. (g). We use SV. — 20 nHz (which is 
the typical value for the Sun) and w = O.OOli? while r^ is 
taken to be the base of the outer convection zone. There 
is no established theory to suggest that the tachocline has 
to coincide with the base of the convection zone, but it 
is quite likely to be true. In some cases, we study the 
variation in the oscillatory signal with the depth of the 
tachocline. The width is chosen to be somewhat small as 
in massive stars the tachocline may be located in regions 
where the scale heights are smaller than that for the solar 
case. The effect of the width of the tachocline on the os- 
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Table 1. Characteristics of stellar models considered 



M/Mq 


Age 


R/Rq 


L/Lq 


Xc 


Vac 


TO 


Convective core 


Convective env. 


He II ionisation 
















rco,o/i? 


Tcorc 


r.„v/i? 


Tcnv 


r^.n/R 


Thoii 




(Gyr) 








(mHz) 


(s) 




(s) 




(s) 




(s) 


0.50 


0.00 


0.460 


0.034 


0.708 


15.092 


1609 






0.632 


1060 


0.977 


320 


0.50 


2.00 


0.479 


0.033 


0.692 


14.093 


1711 






0.614 


1150 


0.976 


350 


0.50 


8.00 


0.484 


0.034 


0.649 


13.759 


1739 






0.611 


1171 


0.976 


355 


0.75 


0.00 


0.640 


0.204 


0.708 


10.218 


2139 


0.115 


2039 


0.687 


1327 


0.981 


396 


0.75 


2.00 


0.666 


0.198 


0.654 


9.560 


2276 






0.670 


1443 


0.980 


429 


0.75 


8.00 


0.689 


0.230 


0.498 


8.837 


2394 






0.664 


1530 


0.980 


459 


1.00 


0.00 


0.861 


0.694 


0.708 


6.929 


2885 


0.147 


2722 


0.724 


1700 


0.982 


516 


1.00 


4.55 


0.976 


0.996 


0.357 


5.321 


3464 






0.721 


2062 


0.980 


645 


1.00 


7.00 


1.071 


1.241 


0.136 


4.399 


3965 






0.715 


2393 


0.979 


763 


1.25 


0.00 


1.159 


1.756 


0.708 


4.585 


4026 


0.140 


3823 


0.789 


2097 


0.983 


684 


1.25 


2.00 


1.325 


2.731 


0.374 


3.417 


4928 


0.053 


4839 


0.847 


2198 


0.983 


797 


1.25 


3.50 


1.531 


3.116 


0.069 


2.616 


6045 


0.053 


5929 


0.791 


3161 


0.979 


1104 


1.50 


0.00 


1.504 


3.777 


0.708 


3.149 


5498 


0.129 


5262 


0.881 


2133 


0.985 


802 


1.50 


1.00 


1.640 


5.786 


0.457 


2.529 


6438 


0.078 


6285 


0.958 


1352 


0.990 


641 


1.50 


1.90 


2.004 


6.340 


0.115 


1.786 


8354 


0.056 


8198 


0.874 


3336 


0.982 


1291 


1.75 


0.00 


1.829 


7.191 


0.708 


2.372 


7061 


0.122 


6794 


0.964 


1367 


0.991 


665 


1.75 


0.90 


2.064 


11.451 


0.331 


1.763 


8566 


0.075 


8375 


0.991 


658 


0.993 


552 


2.00 


0.00 


2.015 


12.505 


0.708 


2.042 


7830 


0.125 


7534 


0.992 


562 


0.994 


485 


2.00 


0.50 


2.064 


19.045 


0.427 


1.720 


8127 


0.091 


7909 


0.991 


628 


0.993 


546 


2.50 


0.00 


2.296 


31.209 


0.708 


1.579 


8609 


0.134 


8257 


0.991 


634 


0.993 


552 


2.50 


0.30 


2.452 


47.999 


0.395 


1.249 


9454 


0.094 


9189 


0.991 


717 


0.992 


623 


3.00 


0.00 


2.546 


65.099 


0.708 


1.301 


9211 


0.143 


8803 


0.992 


637 


0.993 


549 


3.00 


0.20 


2.831 


101.563 


0.364 


0.930 


10707 


0.095 


10401 


0.991 


746 


0.993 


642 


4.00 


0.00 


2.983 


201.651 


0.708 


0.904 


10130 


0.160 


9616 


0.993 


569 


0.995 


475 


4.00 


0.08 


3.079 


295.115 


0.445 


0.825 


10563 


0.119 


10171 


0.994 


561 


0.995 


460 


5.00 


0.00 


3.365 


470.183 


0.708 


0.853 


10873 


0.175 


10255 


0.995 


501 


0.996 


405 


5.00 


0.05 


3.584 


695.338 


0.420 


0.725 


11884 


0.123 


11421 


0.995 


535 


0.996 


415 



cillatory signal has been discussed by Mazumdar & Antia 
(200l|)] For stars with M ^ l.SM© the outer convection 
zone is too shallow and even a width of O.OOli? may be too 
large. But such stars are unlikely to have tachoclines and 
we do not consider this signal in those cases. The varia- 
tion of amplitude in the splitting coefficients with stellar 
mass and age is similar to that in the frequencies. We still 
consider the two separately as the splitting coefficients 
provide a convenient set to study the variation in oscil- 
latory signal with depth, which in this case can be easily 
controlled. 

4. Signal from the base of outer convection zone 

In order to study the oscillatory signature in the frequen- 
cies we calculate the frequencies for each stellar model and 
fit the fourth differences of frequencies to an oscillatory 
signal described in Eq. (0), with two oscillatory compo- 
nents of the same form. The results are shown in Table 0. 
This table also includes the averaged frequency separa- 
tions Ai^o = t'«-(-i,o - I'nfi, Dq = {unfi - i^n-i,2)/6, and 
di/2 = (3/2)(i^„_o — 2i^„_i + i^n+i,o)- The fitted r for the 
two components are close to the actual acoustic depths 
of the He II ionisation zone and the base of the convection 
zone. Thus it is clear that the depths of these layers can be 
estimated from the measured frequencies if the accuracy 
is sufficient to fit the oscillatory signal. 



It can be easily shown that taking the second differ- 
ence of frequencies will enhance the oscillatory signal by 
a factor of 4sin^(27rrAi^), where Az^ is the frequency sep- 
aration {i'n+i,i — Vn/). For t/S.v < 1/12 or rAi/ > 5/12 
this factor will be smaller than unity and the amplitude 
of the oscillatory term will reduce when higher order dif- 
ferences are taken. Thus in such cases it may be better to 
use the second difference or to fit the oscillatory term to 
the frequencies themselves. In the latter case one will need 
some model to fit the smooth part also. It can be easily 
seen that if the smooth part is assumed to be linear in n, 
then it will vanish when second difference is taken. Thus 
we can fit the second difference without worrying about 
the smooth part. Similarly, if the smooth part is assumed 
to be a cubic polynomial in n, it will vanish if the fourth 
difference is taken. For more general smooth functions, 
e.g., splines, the smooth part will not necessarily vanish 
when differences are taken but it will generally become 
smaller. However, the errors also increase as higher order 
differences are computed. Thus, fourth differences may be 
used only if sin2(27rTAi/) > l/^/2 (|Basu et al. 19941). 



If To is the acoustic radius or the sound travel time 
from the centre to the surface of the star, then Ai/ « 
1/(2to). Thus the amplification factor for second differ- 
ence is 4 sin (ttt/tq), which is maximum when r sa tq/2 
or when the layer is halfway between the centre and the 
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Table 2. Results from fits to oscillatory signal 



Frequency separations 



Convective envelope 



Hell ionisation 



M/Mq Age Xc yfp/pQ Auo 



Do 



di 



/2 



Model Sp. coeff. fit Frequency fit Model Frequency fit 



(Gyr) 



(^Hz) (/.iHz) (^Hz) 



(s) 



(s) 



(nHz) 



(s) 



(mHz) 



:to~ 



(s) (s) 



(/iHz) 



0.50 


0.00 


0.708 


2.264 


307.6 


3.124 


18.587 


1060 


1093 


0.114 


1095 


0.096 


320 


315 


0.037 


0.50 


2.00 


0.693 


2.135 


289.6 


2.833 


16.685 


1150 


1183 


0.114 


1185 


0.082 


350 


354 


0.015 


0.50 


8.00 


0.649 


2.102 


284.9 


2.642 


15.580 


1171 


1206 


0.112 


1211 


0.078 


355 


350 


0.033 


0.75 


0.00 


0.708 


1.693 


229.5 


2.926 


17.516 


1327 


1386 


0.112 


1389 


0.088 


396 


386 


0.261 


0.75 


2.00 


0.654 


1.594 


215.8 


2.420 


14.493 


1443 


1505 


0.111 


1509 


0.081 


429 


416 


0.206 


0.75 


8.00 


0.498 


1.514 


205.3 


1.951 


11.463 


1530 


1597 


0.110 


1601 


0.078 


459 


444 


0.203 


1.00 


0.00 


0.708 


1.251 


170.0 


2.512 


15.015 


1700 


1791 


0.117 


1789 


0.069 


516 


531 


0.436 


1.00 


4.55 


0.357 


1.037 


142.1 


1.422 


8.336 


2062 


2176 


0.115 


2181 


0.056 


645 


676 


0.369 


1.00 


7.00 


0.136 


0.902 


124.2 


0.845 


5.948 


2393 


2521 


0.115 


2524 


0.045 


763 


785 


0.301 


1.25 


0.00 


0.708 


0.896 


122.3 


2.066 


12.577 


2097 


2230 


0.121 


2237 


0.047 


684 


747 


0.585 


1.25 


2.00 


0.374 


0.733 


100.4 


1.190 


6.780 


2198 


2315 


0.125 


2336 


0.048 


797 


873 


0.574 


1.25 


3.50 


0.069 


0.590 


81.7 


0.652 


8.481 


3161 


3332 


0.124 


3338 


0.039 


1104 


1114 


0.260 
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Fig. 2. Comparison of amplitude of oscillatory signal in the frequency for evolving stellar models of different masses 
as marked in the figure. All the amplitudes are calculated at a scaled frequency corresponding to 3 mHz for the present 
Sun. a) The amplitude of the signal arising from the base of the convective envelope is shown, b) The amplitude of 
the signal arising from the Hen ionisation region is shown. 



surface in terms of sound travel time. On the other hand, 
if the layer is close to either the surface or the centre, the 
amplification factor will be small. Thus for massive stars 
where the outer convection zone is shallow and the con- 
vective core is also small in terms of acoustic radius, the 
fourth difference will reduce the amplitude significantly. 
In such cases it may be better to use alternate strategy to 
fit oscillatory signal. In this work since we are using exact 
frequencies from stellar models the errors are not ampli- 
fied and hence as long as the amplitude does not decrease 
significantly when differences are taken we can use the 
differences to conveniently suppress the smooth part and 
isolate the oscillatory signal. Even without errors we find 
that when the outer convection zone becomes too shallow 



the oscillatory signal reduces significantly when fourth dif- 
ferences are taken. In such cases we have instead fitted the 
oscillatory signal in second differences. Thus in order to 
facilitate comparison and remove the variation in ampli- 
tude arising from variation of the amplifying factor, we 
reduce all amplitudes to those in frequencies (or splitting 
coefficients) by accounting for the amplification factor in 
taking second or fourth differences. Thus the amplitudes 
given in Table || or those shown in Figs, g, ^ refer to the 
frequencies or splitting coefficients, irrespective of whether 
second or fourth differences were used in the fits. Another 
advantage of scaling the amplitudes to remove the ampli- 
fication factor is that variation in the resulting amplitudes 
is significantly reduced. 
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Fig. 3. Comparison of amplitude of oscillatory signal in 
the splitting coefficients for evolving stellar models of dif- 
ferent masses, arising due to a tachocline at the base of the 
convective envelope. All the amplitudes are calculated at 
a scaled frequency corresponding to 3 mHz for the present 
Sun. Note the difference in y-axis range between this figure 
and Fig. 0. 



Theoretically, there is no reason to believe that taking 
fourth differences of frequencies will remove the smooth 
component. But in practice, using a wide range of stellar 
models, we find that taking the fourth differences essen- 
tially removes the smooth component of the frequency in 
all cases, thus making it much simpler to fit the result- 
ing oscillatory signal. While fitting the second difference 
or the frequencies directly, it is essential to account for 
the smooth part, which necessarily involves some assump- 
tions about how the smooth component is represented. 
The resulting fit will also depend to some extent on these 
assumptions. In general, the smooth part also depends 
on t which makes it more difficult to remove. For more 
rapidly varying signal with r « to/2 like that arising from 
the base of the solar convection zone, it may be relatively 
easy to characterise the smooth part as has been done by 
Monteiro et al. (1994). But when the 'wavelength' of os- 
cillation is comparable to the frequency range in which 
observations are available, it is not easy to distinguish 
between the smooth and oscillatory components. Taking 
the differences also tends to reduce the amplitude in such 
cases making it more difficult to fit the oscillatory signal. 
We believe that this difficulty is more fundamental, as it 
arises from the fact that there is little difference between 
the smooth component and an oscillatory component with 
large wavelength. This problem will manifest itself in some 
form or other irrespective of whether we fit the differences 
or the frequencies directly. If taking the fourth differences 
wipes out the oscillatory signal, then it would be possible 



to approximate it well by a cubic polynomial, which can 
define the smooth part. 

Fig. shows the variation in amplitude of the fitted 
signal in the frequencies with stellar age for a few stellar 
masses. It can be seen that at masses < l.SM© there is not 
much variation in amplitude with age and further even the 
variation with stellar mass is also quite small in the signal 
from the base of the convective envelope. The maximum 
difference being about a factor of two between all models 
included in the figure. The oscillatory signal from the He ii 
ionisation zone has in general larger amplitude and the 
amplitude increases with stellar mass. For M > Mq the 
amplitude of signal from Hen ionisation zone decreases 
with ages, while at lower mass it is essentially indepen- 
dent of age. For more massive stars {M ^ 1.5Mq) the 
acoustic depth of both the He ii ionisation zone and the 
base of convective envelope are comparable and it becomes 
difficult to fit the two signals separately. In addition, be- 
cause the resulting r is quite small, taking the differences 
tend to reduce the amplitudes and it becomes more dif- 
ficult to fit the signal even in the second differences. In 
some cases, when the frequencies of two components are 
close the resulting pattern shows beats in the limited fre- 
quency range that is likely to be available. The apparent 
amplitude of the signal depends on the relative phase of 
the two components and it is difficult to obtain any re- 
liable fit to both the components. The characteristics of 
the signal arising from the base of the convective envelope 
may also be different in these stars as the outer convec- 
tion zone splits into two parts. Consequently, wc have not 
shown the results for these stars. 

Apart from the frequencies we also calculate the split- 
ting coefficients for each of the stellar models with a model 
rotation profile given by Eq. (||). This rotation profile as- 
sumes a tachocline at the base of the outer convection 
zone. The fourth difference of the calculated splitting co- 
efficient C3{n,£) for £ = 2,3 is again fitted to an oscilla- 
tory signature of the form given in Eq. (|j) to calculate the 
characteristics of the tachocline. These results are also in- 
cluded in Table g and it is clear that the fitted r is close 
to the actual depth of the tachocline in the model profile. 
The variation of amplitude with stellar mass and age is 
shown in Fig. ^. In this case, the variation with stellar 
mass and age is very small. 

The typical amplitudes of the oscillatory signal in the 
frequencies due to the base of outer convection zone is 
0.05 /iHz, while that due to He ii ionisation zone is 0.5 /xHz. 
We have not added any errors in the model frequencies 
and hence there is little difficulty in identifying this sig- 
nal for low mass (M < 1.5Mq) stars. Actual observa- 
tions will naturally have some random error associated 
with each frequency. From simulations with artificial data 
sets it turns out that it is possible to detect the oscillatory 
signal as long as the errors in frequencies are less than or 
of the order of the amplitude of oscillations. Thus we need 
an accuracy of about 0.05 ^Hz to detect this signal from 
the base of convective envelope and measure its character- 
istics. For T « To/2 the signal to noise ratio will improve 
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by a factor of two when fourth differences are taken and 
an accuracy of about 0.1 /iHz may be sufficient to measure 
the characteristics of the signal. Such an accuracy may be 
possible with the forthcoming space missions like COROT 
(iBaglin et al. 1998D, MOST ([Matthews 1998|) and MONS 



(Kjeldsen & Bedding 1998). The oscillatory signal due to 
the He ii ionisation zone has an order of magnitude larger 
amplitude for solar mass stars and it would be easier to 
detect this signal. This signal may be useful in measuring 
the helium abundance in stellar envelopes. The oscillatory 
signal in splitting coefficients has much smaller amplitude 
and unless the stars are rotating much faster and have a 
much larger S^ it may not be possible to detect such a 
signal. 



5. Signal from the edge of convective core 

5.1. Signal in frequencies of iiigii mass stars 

For stars with AI <; 1.2M0, apart from the outer con- 
vection zone, a part of the inner core is also convective 
and we would expect an additional oscillatory component 
in the frequencies arising from discontinuity in derivatives 
of sound speed at the outer boundary of the convective 
core. In this case, t fn tq and as mentioned in the previ- 
ous section, taking fourth difference will reduce the ampli- 
tude significantly. In order to understand the oscillatory 
signal, we consider a stellar model for M = 2Mq at an 
age of 5 X 10* years. This model has a convective core 
with acoustic radius COStq and taking the fourth differ- 
ence will reduce the amplitude by a factor of about 1000. 
The outer convection zone in this model has an acous- 
tic depth of about O.OStq and that oscillatory signal will 
also reduce by a factor of 20. Fig. Q shows the second and 
fourth differences of frequencies for this model. It is clear 
from the figure that the oscillatory signal arising from the 
outer boundary of the convective core is not seen in either 
the second or the fourth differences. These differences do 
show the signal from the base of the outer convective zone 
or the Hen ionisation region, while the oscillations that 
are seen in the fourth differences have completely differ- 
ent T. Fitting an oscillatory signal in a restricted frequency 
range shows that this signal is arising from r ~ 5300 s, 
which corresponds to r/ R ^ 0.6, where we do not expect 
any boundary of convective region. 

In order to understand this signal we show in Fig. ra, 
the quantity W{r) — {1 / g)dc'^ / dr for a few different stellar 
models. It can be seen that while low mass stars do not 
have any sharp variation in this scaled derivative of the 
sound speed in the radiative interior, the high mass stars 
all show a prominent dip in W(r) in the radiative interior. 
This dip is quite small compared to the sharp rise near the 
base of the outer convection zone in near-solar mass stars. 
But for high mass stars where the outer convection zone is 
very shallow, this dip becomes the dominant feature in the 
first derivative of the sound speed and hence contributes to 
the oscillatory signature. Moreover, calculating the fourth 
difference increases the amplitude of this signal by more 
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Fig. 4. Second (upper panel) and fourth [lower panel) dif- 
ferences of frequencies for a 2Mq star at an age of 0.5 Gyr. 
Only modes with £ = to 3 have been considered. Fits to 
the data with a function of the form described in Eq. ^ 
have been shown for the full range of frequencies in case of 
the second differences, and for a restricted range in case of 
the fourth difference. Note that the range of y-axis in this 
figure is an order of magnitude less than that in Fig. n. 

than an order of magnitude, while amplitudes of other 
oscillatory terms coming from edges of convective layers 
are significantly reduced. Thus in the fourth difference the 
prominent oscillations are due to this dip in W{r), while in 
the second difference the oscillations arising from the base 
of the outer convection zone are seen. But the expected 
signal from the outer edge of the convective core is not seen 
in either the second or the fourth differences. Since there 
is no discontinuity in the derivatives of the sound speed in 
the radiative region the oscillatory signal does not have the 
same form as that expected from a discontinuity, but for 
low frequency modes, where the radial wavelength would 
be larger than the size of the dip the signal is seen more 
clearly. At high frequencies the radial wavelength may be 
comparable to or smaller than the size of the dip in W{r) 
and the oscillatory signal is not clear. Further, just like 
the He ii ionisation zone, the amplitude of this oscillatory 
signal should also be modulated by an oscillatory factor, 
which is not included in our fits. 

5.2. Signal in splitting coefficients 

In order to understand our failure to see any oscillatory 
signal from the convective core in the frequencies, we at- 
tempt to study the oscillatory signal due to the tachocline 
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Fig. 5. The quantity W — {l/g)dc^/dr for stars of dif- 
ferent masses and ages. 

Table 3. The variation of amplitude of S'^c^ with position 
of tachocline for a 2Mq star of age 0.5 Gyr. All the am- 
plitudes are calculated at the scaled frequency 1.19 mHz, 
corresponding to 2.5 mHz for a solar model. The total 
acoustic depth of the star is tq — 8128 s. 
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in the splitting coefficients. We calculate the splitting co- 
efficients assuming the tachocline to have different depths 
to study the variation in amplitude with depth and the 
results are shown in Table |^. We have considered a stellar 
model with mass 2Mq and age 0.5 Gyr for the purpose of 
illustration. The amplitude appears to be maximum when 
the tachocline is located near the acoustic midpoint of the 
star (in this particular case, near r/R = 0.8) and falls off 
as it shifts to either end. It can be easily seen that this 
variation is more or less fully explained by the amplifica- 
tion factor while taking the differences. Thus it appears 
that the amplitude of the signal in the splitting coefficients 
is more or less independent of the depth and the apparent 
variations are coming from calculating the differences. It is 
clear that in order to detect the oscillatory signal from the 
outer edge of the convective core, we should look at the 
splitting coefficients directly, rather than their differences. 

Fig. g shows the second difference of splitting coeffi- 
cients for ^ = 2 modes arising from the tachocline at dif- 



ferent depths in the stellar model mentioned above. It can 
be seen that the wavelength of oscillations reduces starting 
from rd = 0.95R as we go deeper as expected. However, 
after r^ — 0.8R the wavelength appears to increase again 
contrary to what we expect. This is similar to the phe- 
nomenon of aliasing in Fourier transform. For simplicity 
if we assume that the frequencies of £ = 2 modes are uni- 
formly spaced then the oscillatory term can be written 
as 

C3{n,i) = Asin[4:TTT{i^o + nAi'f) + cf)] 

— ylsin[47rTi'o + + 27rnr/To], (5) 

where Ai>i is the spacing between successive modes for 
given £. Here we have used the fact that for low £, 
Ai^i « 1/(2to). Now if t = tq it is clear that for all n the 
oscillatory term will be the same and C3(n, i) will not show 
any oscillations. If we take t' = tq — t, then it is clear that 
the oscillatory term will be the same apart from a change 
in the phase for all modes. Thus from the oscillatory sig- 
nal for an individual £ it is not possible to distinguish 
between these two values of r. Fig. ^ shows the actual fit 
obtained using two different r for some depths. This is 
exactly similar to aliasing in discrete Fourier transform. 
Alternatively, it means that we can measure the acoustic 
distance from either boundary. If we consider sampling at 
a uniform spacing of Av « 1/(2to), then the correspond- 
ing Nyquist frequency will be l/{2Aiy) « tq. Thus 'fre- 
quencies' higher than tq will be aliased to lower frequen- 
cies. Because of the factor of 2 inside the sine function in 
Eq. (0) this actually corresponds to an acoustic depth of 
To/2. Thus this is identical to aliasing in discrete Fourier 
transform. 

It may be noted that the lowest panel in Fig. ^ for 
Td = O.IR is shown separately, as in this case instead of 
the second differences (which have very small amplitude), 
we show the splitting coefficients themselves. Although, 
the amplitude of variation in this is comparable to am- 
plitudes at other depths, it is impossible to distinguish 
this signal from a smooth trend. At other depths we have 
shown the second difference as the presence of smooth 
trend in splitting coefficients does not show the oscilla- 
tory signal clearly. It is clear from this panel that both 
i = 2 and 3 modes separately fall on a smooth curve and 
only if they are combined together the oscillatory signal 
from the core can be identified. It may be noted that the 
form of rotation rate used to calculate these splitting co- 
efficients does not have any smooth component and hence 
the smooth component in resulting splitting coefficients 
is quite small, thus enabling us to see the signal in some 
form. In real stars there will also be some smooth compo- 
nent which will be difficult to isolate from the oscillatory 
part when the apparent r is quite small. 

In principle, the ambiguity in r can be resolved by 
including multiple £ values. The two panels at the lower 
end of Fig. |g shows the splitting coefficients for r^ — 0.3R 
and O.li? for £ ~ 2,3. It is clear that although modes 
corresponding to each £ show a variation with large wave- 
length, the two do not agree with each other in phase 



Mazumdar and Antia: Seismic study of stellar convective cores 



N 







N 

X 

c 



1.0 1.2 

!/ (mHz) 

Fig. 6. Second difference of the splitting coefficient C3 for 
£ = 2 only {filled circles), for a 2Mq star at an age of 
0.5 Gyr, arising out of a tachocline at different depths, 
rd as marked in each panel. A fit to the data using a 
function of the form described in Eq. (jj) is shown in each 
case. The two panels at the bottom shows the fit with the 
"true" value of r^ (dotted line), which is obtained when 
both £ — 2 and 3 [empty circles) are considered. Compare 
this to a fit with the "aliased" Td (solid line), obtained 
with £ — 2 modes only. 



and if we attempt a fit using Eq. (Q), the fit will not be 
good when r corresponds to the smaller value implied by 
individual £. Only an oscillatory form with larger value 
of T can fit all the modes. In terms of discrete Fourier 
transform this can be understood as follows: because of 
additional £ the sampling interval is reduced thus increas- 
ing the Nyquist frequency, enabling us to determine higher 
frequencies. However, the problem with fitting any form 
is that in all cases there will be some smooth component 
in the splitting coefficient arising from smooth variations 
in the rotation rate and it will be difficult to separate out 
this component which has to be done for each £ separately, 
as the smooth component will depend on £ and frequency. 
The problem arises because in this case the oscillatory 



component also appears to be smooth due to aliasing and 
it will be difficult to isolate it from the smooth component 
as can be seen by the panel for r^ = OAR in Fig. a This 
separation will be more difficult when errors are added to 
the splitting coefficients, as will be the case in real ob- 
served values. Thus it is difficult to detect any possible 
oscillatory signature arising from convective cores in high 
mass stars. 



5.3. Signal from the convective core 

Detecting the oscillatory signal in the frequency from the 
convective core is even more difficult as the frequencies 
have dominant oscillatory contributions arising from other 
variations in sound speed derivatives. Thus the signal due 
to the outer convection zone will have similar r to that 
from convective core (after aliasing) and it will be difficult 
to isolate the core contribution from the surface contribu- 
tion. Further, there is a smooth variation with n for the 
frequencies, which is several orders of magnitude larger 
than the oscillatory part and it will be more difficult to 
remove this smooth variation. Of all the models tried we 
could detect this signal in only 1.75Mq star with an age 
around 0.9 Gyr (Fig. |^). Interestingly even with exact fre- 
quency data this signal is not clearly seen at all ages. We 
believe the main reason for this is that at this particular 
age the amplitude of signal arising from the outer con- 
vection zone and Hell ionisation zone combine to yield 
a small amplitude, possibly because of beating between 
the two frequencies which are reasonably close. Because 
of the small amplitude, even smaller oscillatory signal due 
to the convective core shows up with amplitude of order of 
0.01 /iHz, which would otherwise have been missed. This 
signal will translate to an amplitude of 0.3 /iHz in frequen- 
cies, which can be compared with frequency separation of 
50 /LtHz. Although, the amplitude of this oscillatory signal 
is fairly large in the frequencies themselves it gets reduced 
when differences are taken. If we do not take differences 
we need to remove the smooth part independently, which 
is also difficult as explained earlier. We also tried to use 
the first difference of frequency instead of second differ- 
ence, so that the amplitude of oscillatory signal from the 
convective core would be larger. But even in that case the 
smooth part is substantial and it is not easy to remove 
it in most cases, thus making it difficult to identify the 
oscillatory signal. 

In an attempt to detect the oscillatory signature we 
consider even higher mass stars where the convective core 
is comparatively larger. Thus we consider a stellar model 
with M = 5Mq and age 5 x 10^ years. In this case the 
convective core has a radius of 0.04ro and Fig. M shows 
the second differences of the frequencies. It appears that 
at a somewhat high frequency range some oscillatory sig- 
nal which is similar to what one would expect from the 
convective core is present. However, the acoustic cutoff 
frequency for this model is only 0.72 mHz and hence the 
oscillatory pattern is seen only in high frequency modes 
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Fig. 7. Second difference of the frequencies, for a 1.75Mq 
star at an age of 0.9 Gyr. A fit to the data using a function 
having two similar components as described in Eq. (jj) is 
shown. The high frequency oscillations are due to outer 
boundary of the convective core, while the dominant os- 
cillation is due to a combination of Hen ionisation zone 
and the base of the convective envelope. 

which are not trapped in the interior. These modes are un- 
likely to be seen in distant stars. Even if they are seen, it 
will be difficult to determine their frequencies accurately, 
because of large widths in the power spectrum that may 
be expected. Similar signal at frequencies above the acous- 
tic cutoff has been seen in some other stellar models also. 
The amplitude of dominant oscillations due to the Hen 
ionisation zone falls off rapidly with frequency and that 
probably allows the smaller signal arising from the con- 
vective core to be seen at high frequencies. 

Since we have failed to detect the oscillatory signal 
from convective core, in the next section we propose an 
alternative technique which looks at the variations in the 
smooth part to detect convective cores in stars and possi- 
bly measure its size. 

6. Signature of convective core in frequency 
separation 

From the asymptotic theory of stellar oscillations it is well 
known that the frequency separations 






Dp, = 



Al + Q 



(6) 



^^1/2 — :^(^n,0 ~ 2i^„.i + fn+i^o). 



are related to variation in the sound speed across the stel- 
lar radius including the convective core. Some of these fre- 
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Fig. 8. Second difference of the frequencies, for a ^Mq 
star at an age of 0.05 Gyr. A fit to the data in a restricted 
frequency range using a function having two similar com- 
ponents as described in Eq. (0) is shown. 
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Fig. 9. Evolutionary tracks for different stellar masses in 
the (AfQ, I?o) diagram. Both the frequency separations 
have been scaled by ^JpJpQ. 

quency separations are quite sensitive to conditions in stel- 
lar core. These frequency separations can be easily mea- 
sured for stars using only the low degree modes. Here we 
consider only the frequency separations for i? = aver- 
aged over a frequency range scaled appropriately for each 
model to correspond to a range of 3-4 mHz for the Sun. 

Ghristcnsen-Dalsgaard (1993) has suggested that these 
frequency separations can be related to the stellar age and 
mass from the evolutionary tracks in the (Ai/Qj-Do) plane. 
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Fig. 10. Large frequency separation, Avq, as a function 
of convective core radius. Both the axes have been scaled 

by \JpIpq- 



which is also known as the C-D diagram. Since the fre- 
quencies of stellar oscillations are expected to scale as 
yjMjB? it will be more instructive to scale these fre- 
quency differences by this factor to isolate other effects 
in the frequency variation. Fig. ^ shows the scaled fre- 
quency separations plotted against each other for vari- 
ous stellar models. This is similar to the figure shown by 



Christenscn-Dalsgaard (1998). It can be immediately seen 



that after scaling the variation in the frequency separa- 
tion is significantly reduced and moreover, all models for 
M ^ 1.25M0 appear to have the scaled large frequency 
separation Avo/y/p/pQ « 137 ^Hz, close to the solar 
value. These are the models that do not have any sig- 
nificant convective core. All models with convective core 
appear to have significantly smaller scaled frequency sepa- 
rations. Thus it appears that if the stellar mass and radius 
are known then by looking at the scaled frequency sepa- 
ration we can infer the presence of a convective core (cf., 
Thompson 200C| ) . If the scaled large frequency separation 
is less than about 134 /iHz, then the star is likely to have 
a convective core. 

Apart from detecting the presence of convective core 
we would also like to measure its size. Indeed, the asymp- 
totic theory of stellar oscillations (e.g., Roxburgh & 



Vorontsov 2001) may be able to provide a relation con- 
necting the frequency separation to the size of the con- 
vective core. In this work, using the frequencies computed 
from different stellar models we attempt to check if the 
deviation from the solar value of the large frequency sepa- 
ration can be correlated to the size of the convective core. 
It might also be possible to derive some combination of the 
large and small frequency separations which would vary in 
tandem with the size of the convective core. In Fig. M we 



have plotted the difference of the scaled large frequency 
separation, from the solar value of 137 /iHz, against the 
radius of the convective core for stellar models of different 
mass and age. Both the axes have been scaled to bring 
the curves corresponding to different masses close to each 
other. We find that there is an approximately linear re- 
lation between these two quantities for a wide range of 
masses. Although, all the curves do not collapse into one 
single relationship, this diagram might be used to esti- 
mate the size of the convective core, provided the mass 
and radius of the star are independently known. Only the 
curves for 1.25Mq and l.SM© stars appear to fall outside 
the band. 

The small frequency separation, Dq, is found to bear 
a correlation with the radial size of the convective core. 
This relation becomes clear with appropriate scaling of 
Dq as well as Tcoic, as shown in Fig. [ll](a). One problem 
with using Dq to measure the convective core size is that 
the range of values covered by a solar mass star during its 
evolution has a large overlap with the value for stars with 
convective core. Thus by looking at the value of Dq alone 
it will be difficult to determine the convective core size. 

Basu et al. (2001) have pointed out that the behaviour 
of the small separation, di/2, is somewhat different from 
the other frequency separations. We find that this quan- 
tity does not have a clear relationship with the size of the 
core. Nevertheless, as shown in Fig. 0(b), its value derived 
for models with a convective core (after some scaling) is 
conspicuously different from that of models without con- 
vective core. However, for both the small separations, we 
are unable to collapse all the curves together as would 
have been desirable for them to be used as definitive di- 
agnostics for the size of a convective core. However, they 
can be used as complementary diagnostics for properties 
of convective core. 

An advantage of using the frequency separation as op- 
posed to the oscillatory signal discussed in the previous 
section is that comparatively lower accuracy is required 
in this case. Thus an accuracy of order of 1 /iHz may be 
sufficient to get an estimate of convective core size using 
the frequency separations. On the other hand, extract- 
ing the oscillatory signal will require at least an order 
of magnitude higher accuracy in measured frequencies. 
Even with very high accuracy it will be difficult to iso- 
late the oscillatory signal from convective core from other 
signals. However, the use of frequency separation requires 
the knowledge of stellar mass and radius for scaling the 
frequencies. This may not be a problem for nearby stars 
which are likely to be the prime targets of asteroseismic 
missions. 

7. Conclusions 

In this work we have tried to study the oscillatory signal 
in the frequency arising from discontinuities in derivatives 
of sound speed inside a star. This oscillatory signal can 
be conveniently used to study the location of the base of 
the outer convection zone or the Heii ionisation zone as 



12 



Mazumdar and Antia: Seismic study of stellar convective cores 




0.00 0.02 0.04 

(p/PG)(^core/«) 



0.06 



3. 



o 

IQ. 
\ 







0.00 0.02 0.04 

(p/Po){^core/«) 



0.06 
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has been suggested by Monteiro et al. (200C ) . Our results 
also support their conclusions. We have studied this sig- 
nal in a wide range of stellar models covering a range of 
stellar masses and ages. One of the reasons for this study 
was to possibly identify promising stars where such signal 
may be large enough to be detected. However, from our 
study it appears that there is only a small variation in the 
amplitude of the oscillatory signal arising from the base 
of the convection zone between stars of different masses 
and ages on the main sequence. The amplitude increases 
with decreasing mass or age. The amplitude of the oscilla- 
tory signal arising from the He ii ionisation zone appears 
to increase with stellar mass and can possibly be used to 
measure helium abundance in stellar envelopes. Typical 
amplitude of oscillations from the base of the outer con- 
vection zone is around 0.05 /iHz and from the simulations 
that we have carried out it appears that comparable accu- 
racy will be required to extract this signal and measure the 
characteristics of the convection zone boundary. Although, 
in this study, we have not included any overshoot below 
the convection zone, it is well known that the amplitude of 
this signal increases with overshoot and it will be possible 
to use this signal to measure the extent of overshoot b elow 
the base of the convection zone ( Monteiro et al. 2000 ). For 
the case of the Sun, this has already been done, though 
with the help of higher degree modes whose frequencies 
are known more accurately. The oscillatory signal arising 
from the He ii ionisation zone is an order of magnitude 
larger and it should be possible to detect this signal in 
observed frequencies. 

In higher mass stars we find an oscillatory signal aris- 
ing from the dip in Wir) in the radiative interior, which 
may also provide some diagnostic for stellar models. From 
Fig. 1^ it can be seen that the location of the dip depends 



on stellar mass. Since this location can be determined 
through T in oscillatory signal it can provide some con- 
straints on stellar models. It can be shown that in radiative 
regions W(r) is also affected by opacity, which in turn will 
depend on heavy element abundance, Z . Thus the depth 
of the dip also depends on Z . However, the amplitude 
of oscillations arising from this dip is about 0.05 /iHz in 
the fourth differences and hence we will need accuracy of 
order of 0.01 /iHz in measured frequencies to detect this 
signal. This level of accuracy may not be achieved in im- 
mediate future, but with further improvements and longer 
observations it may be possible to detect this signal. 

Most of the earlier studies have concentrated on the os- 
cillatory signal from the outer convection zone. It is clear 
that there would be a discontinuity in the derivative of 
sound speed at the outer edge of the convective core in a 
massive star. This discontinuity should also give rise to a 
similar oscillatory signal. Although, Monteiro et al. (1998 ) 
have suggested that this signal can be used to study char- 
acteristics of convective core, it is not clear if this signal 
has been seen in frequencies of a stellar model. Thus in 
this work we attempt to study the oscillatory signature 
arising from a convective core. Even using exact frequen- 
cies from stellar models we find it difficult to detect this 
signal. We have tried to identify the reasons for this fail- 
ure by studying a similar signal in the splitting coefficients 
arising from a possible tachocline at different depths. We 
find that the oscillations display a phenomenon similar to 
aliasing in discrete Fourier transform, because of which 
using modes from a single value of ^, it is not possible to 
distinguish between oscillations arising from discontinu- 
ities at acoustic depths of r and rg — r, where tq is the 
total acoustic radius of the star. Due to this aliasing the 
signal from the convective core would appear as a smooth 
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variation in frequencies from any individual L Though, in 
principle, by combining data for ^ = 0,1,2,3 we should be 
able to characterise this signal, but in practice we have to 
remove the smooth part for each I before the oscillatory 
part can be studied. It is not easy to separate this signal 
from the convective core, thus making it difhcult to study 
its characteristics. 

From our attempts to find the oscillatory signal from 
the convective core, we could find this signature in a few 
stellar models at frequencies above the acoustic cutoff fre- 
quencies. Such modes are not likely to be observed in 
stars, but from these fits it appears that the amplitude of 
the oscillatory signal in frequencies is actually quite large, 
around 0.3 /iHz, but because of our inability to separate 
it from the smooth component it is difficult to detect this 
signal. Of all the models that we tried we could fit this 
signal only in a 1.75M0 star around the age of 0.9 Gyr at 
frequencies less than the acoustic cutoff frequency. This is 
most probably because the dominant oscillatory signal due 
to outer convection zone and He ii ionisation zone combine 
to yield a small amplitude in this stellar model, because 
of beating. Thus we conclude that an oscillatory signal of 
the form studied in this work does not provide a viable 
means of studying convective cores in massive stars. 

Having failed to detect the oscillatory signature from 
the convective core in frequencies of stellar oscillations, 
we attempt to check if the smooth part of the frequency 
can be used to study the properties of the core. The fre- 
quency separations Ai/^, D( and dii2 defined by Eq. (R) 
are known to be sensitive to stellar structure including 
the core (e.g., phristensen-Dalsgaard 199S ; Basu et aL| 
200l[)T]Since the frequencies of stellar oscillations gener- 
ally scale with ^MjB? it is instructive to consider the 
scaled frequency separations. It turns out that the scaled 
large frequency separation is generally around 137 /iHz for 
stars without convective core, while for massive stars with 
convective cores it tends to be significantly smaller. Thus 
this scaled frequency separation can be used to detect the 
presence of a convective core. After suitable scaling the de- 
parture of this frequency separation from the mean value 
for stars without a convective core can be roughly cor- 
related to the size of the convective core as can be seen 
from Fig. |lO|. It can be seen that apart from intermediate 
masses of 1.25M0 and 1.5M0, for higher masses all the 
evolutionary tracks lie in a band and from the measure 
of the frequency separation together with the scaling fac- 
tor it should be possible to estimate an approximate size 
of the convective core. The small frequency separations 
I?o and d\i2 are also sensitive to properties of the core, 
though we could not find a similar tight relation with the 
core size. We believe that by combining all the informa- 
tion in the smooth component of the frequency for only 
low degree modes it should be possible to study some prop- 
erties of the convective core, including its size. An advan- 
tage of using these frequency separations is that compara- 
tively lower accuracy of 1 /iHz in frequency would be suffi- 
cient, while we would need much higher accuracy of about 
0.05 /xHz to study the oscillatory signal. Even the oscil- 



latory signal from the base of the solar convection zone 
has been mostly studied using intermediate degree modes 
which have higher accuracy. In this work we have consid- 
ered only the mean frequency separations which naturally 
discards some information, which is contained in the indi- 
vidual frequency separations for each value of n. In partic- 
ular, the variation of frequency separation with n will also 
give some useful information about the convective core. 
Instead of frequency separation inversion techniques us- 
ing only low degree modes have also been tried to study 



stellar cores (Basu et al. 2001). 

In the foregoing discussion it is assumed that the scal- 
ing factor ^ MjB? is known for stars. This would require 
knowledge of stellar mass and radius. In principle, the stel- 
l ar mass and age can be deter mined from the C-D diagram 
( Christensen-Dalsgaard 1993 ), while if the distance to the 
star is known, its luminosity and hence the radius can be 
determined. From Fig. ^ it can be seen that the scaled 
large frequency separation is not very sensitive to stellar 
mass or age. The variation over the entire range of stel- 
lar models considered by us is less than 10%. Hence in 
absence of any other information the measured large fre- 
quency separation can also be used to measure the scale 
factor. But in that case the resulting separation cannot be 
used to measure the size of convective core. 
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